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ABSTRACT 


This  report  presents  an  esq^erinental  study,  throu^  pressure  dlstrlb* 
utions,  of  the  honogeneous  nucleation  and  condensation  of  initially  super¬ 
heated  aauBonia  vapor  during  its  expansicm  in  a  low  Mach  nuaber  supersonic 
nozzle.  The  results  of  the  e3q>eriiBent8  are  ccs^ared  with  the  classical 
nucleation  theoxy  of  Becker  and  Coring  or  Frenkel-Zeldovich,  on  the  basis 
of  the  predicted  pressiure  and  temperature  at  which  condensation  would  occur 
starting  from  a  given  stagnation  condition. 

It  is  found  that  ammonia  shows  far  less  supersaturation  than  this 
theoxy  predicts.  The  effect  of  arbitrary  corrections  to  the  surface  ten¬ 
sion  showed  that  a  further  reduction  in  surface  tension  of  the  drop  relative 
to  the  flat<<fila  value  gives  better  agreement  between  measured  and  predicted 
pressiure  distributions .  Even  better  agreement  is  obtained  without  a  sur¬ 
face  tension  correction  if  we  assume  that  the  nucleation  rate  is  greater 
17 

by  a  factor  of  10  .  This  seems  to  lend  some  support  to  the  recent  theories 

of  Lothe  and  Pound^^^^and  Oriani  and  Sundquist^^^ ^concerning  gassification 
corrections  to  the  free  energy  of  formation  of  critical  droplets.  It  is 
not  certain  however  how  much  of  this  correction  can  be  assigned  to  uncer¬ 
tainties  in  surface  tension,  and  how  much  to  the  'gassification'  effect. 
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CONDENSATION  OF  AMiONIA  VAPOR  DURING  RAPID  EXPANSION 


by 

Michael  Kreauner  &  Olufemi  Oktirouxunu 

I.  INTRODUCTION 

A»  Background  to  the  Problem 

Interest  in  the  phenonenon  of  condensation  of  vapors  undergoing  various 
forms  of  ej^ansion  have  gained  considerable  momentum  amongst  scientists  and 
engineers  over  the  last  fev  decades.  HQvever,  the  state  of  knovledge  in  this 
field  over  the  same  period  has  in^roved,  if  at  all,  by  ▼exy  little  and  is  at 
the  moment  in  a  rather  imsatlsfactory  condition.  !lhis  is  not  to  say  that 
no  progress  has  been  made,  indeed,  several  heretofore  unrecognized  factors 
that  do  have  a  significant  Influence  on  condensation  have  nov  been  emphasized 
by  recent  Investigators,  but  to  the  extent  that  there  is  not  in  general  a 
consensus  of  expert  opinion  on  the  exact  nature  of  the  interrelationships 
between  these  variables,  the  role  of  experiments  cannot  be  over'-emphasized, 
more  or  less  as  an  arbiter  between  the  various  existing  theories  of  nucle- 
ation. 

By  condensation,  in  this  report,  we  refer  to  that  resulting  from  homo¬ 
geneous  nucleation  from  the  vapor  phase,  the  nuclei  in  this  case  being  clus¬ 
ters  of  atoms  of  the  vapor  which  increase  in  number  as  the  supersaturation 
of  the  vapor  increases,  and  subsequently  cause  a  collapse  of  the  super¬ 
saturated  state . 

Some  of  av  a'^liest  xperiments  on  this  phenomenon  were  conducted  by 
Wilson^^^.  ttud  later  by  Volmer  and  71ood^^\  all  utilizing  cloud 

chambers  for  the  expansion,  although  a  fev  earlier  nozzle  experiments  had 

(27) 

been  done  by  Him  and  Cazln^  '  .  It  is  well  Known  however  that  a  cloud 
chamber  expansion  is  inherently  an  unsteady  process,  is  very  much  influenced 
by  the  presence  of  dust  or  any  other  foreign  particles,  since  the  expansion 


2 


rate  is  very  sloe  and  hence  allows  considerable  interaction  between  the  dust 

particles  and  the  vapor.  Besides,  observation  of  condensation  in  cloud 

chaaber  eiqperiaents  depends  on  visual  observation  of  the- appearance  of  a 

(1*) 

cloud,  idtereas  Hill  ,  by  coiqparing  Yellot's  separate  neastirenents  of  vis* 
ible  cloud  and  pressure  distributions,  has  concluded  that  the  actual  appear¬ 
ance  of  visible  clouds  occurs  a  little  later  than  the  beginning  of  conden¬ 
sation.  Since  the  entire  process  is  unsteady,  it  is  impossible  to  obtain 
a  history  of  the  events  prior  to  condensation.  One  could,  therefore,  wonder 
whether  cloud  chaaber  data  can  truly  be  used  to  check  existing  nucleati(m 
theories,  mostly  based  on  the  assunption  of  c  steady  state  distribution  of 
nuclei;  and  one  might  say  that  any  apparent  agreement  between  theory  and 
data  obtained  from  cloud  chambers  would  at  least  be  partly  fortuitous  and 
would  not  truly  reflect  the  validity  of  the  theory  itself. 

Other  experimenters  have  studied  the  same  phenomenon  through  nozzle 
expansions,  using  superheated  steam  as  the  vapor.  The  data  of  Stodola^^\ 
Xellot^^^*^*^^  fietalllata^®\  Binnie  and  Woods^^^  and  those  of  Wegener  & 
Pouring^ are  well-known  in  the  literatxure.  Hozzle  expansions  are  no  doubt 
more  reliable  than  cloud  chamber  sxpa.'.sions ;  steady  state  cim  very  readily 
be  maintained,  and  because  of  the  very  high  expansion  rates,  dust  paxiticles 
cannot  have  any  appireciable  effect  on  nucleation.  It  has  been  shown  ly  Oswa- 
titsch^^^^that  dust  pax^icles  with  a  concentration  of  as  high  as  10^  part¬ 
icles  per  cc.  would  have  virtually  no  effect  on  the  condensation  of  steam 
in  a  nozzle.  The  nozzle  data  mentioned  above  conpairr's  favourably  with  the 

(li) 

Becker-Doring  equation,  as  shown  by  Hill  .  But,  since  all  the  available 
data  in  the  literature  is  on  steam,  it  would  seem  that  a  true  test  of  excist- 
ing  theories  cannot  be  made  until  reliable  data  on  other  vapors  becomes 
available.  It  is  with  this  in  mind  that  the  present  experiments  have  been 
conducted.  It  is  hoped  thnf  it  v/ill  supplement  excisting  data  on  steam  in 
enabling  us  to  determine  Just  how  useful  existing  theories  are  in  1)  pre- 
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dieting  the  onset  of  condensation,  and  2)  describing  the  histoxy  of  the  con*' 
densation,  particularly  in  terns  of  pressure  distributions  in  the  vapor, 
both  before  and  after  the  aettial  condensation  process. 

It  will  not  be  out  of  place  at  this  point  to  ask  Just  vhy  there  is  a  need 
for  a  better  understandinf^  of  such  a- -single  process  as  condentatiOn  im^a  .vapor? 
While  the  physicist  or  chenist  may  be  Interested  for  qxiite  different  reasons, 
the  engineer  is  primarily  concerned  vith  the  effects  of  possible  condensation 
of  the  working  fluid  in  a  power  plant  on  the  cycle  efficiency  and  other 
uesign  parameters.  Since,  in  the  Rankine  cycle  with  saturated  vapor  at  the 
turbine  inlet,  condensation  of  the  working  flxxid  is  almost  inevitable,  a 
good  design  would  have  to  take  this  into  c<xisideration,  by  making  allowance 
foi'  losses  due  to  a)  the  drag  of  the  condensed  droplets  and  b)  dre^  inpinge- 
ment  on  successive  rows  of  blades .  The  exit  are;i  required  to  pass  a  given 
flow  would  also  have  to  be  greater  if  condensation  occurs  in  the  working 
fluid. 

B.  Hucleation  Itoeoiy 

There  have  been  several  forms  of  the  nucleation  rate  aquation  developed 

(12) 

by  various  authors,  outstanding  among  whom  are  Volmer  and  Weber  ,  Becker 

(13)  (ll*) 

and  Boring  and  Zeldovlch  .  An  outline  of  the  derivation,  together 
with  the  necessary  assumptions  is  presented  below,  along  the  lines  of  Frenkel, 
Becker,  and  Boring 

It  is  assumed  that: 

1)  there  are  present  in  the  Vflq)or,  at  any  instant  of  time,  micro¬ 
scopic  embryos  of  the  liquid  phase,  with  varying  sizes  acting  as  nuclei, 
and  differing  fr<xtt  the  macroscopic  state  of  the  liquid  only  in  their 
size. 

2)  The  thermodynamic  state  of  these  embryos  may  be  described  by 
macroscopic  properties  of  the  liquid  phase,  even  though  they  contain, 
in  general,  only  a  few  molecules. 


u. 

3)  The  energy  of  fonuitlon  of  the  enbxyo  consists  of  two  suns: 
a)  the  Gibbs  free  energy  change  associated  with  the  condensati<»  of 
bulk  liquid  froa  the  vapor  phase,  and  b)  the  surface  energy  associated 
with  the  coaainution  of  the  bulk  liquid  into  spherical  droplets. 

Thus,  the  free  energy  change  in  the  formation  of  a  g-sized  drqp  is 
given  hy 

AI  =  (*!*•-  +  4irrV 

idiere  ■  Gibbs  free  energy /molecule  in  the  liquid  phase 

•  Gibbs  free  energy /molecule  in  the  gaseous  phase 
^  «  number  of  molecules  in  the  drqp 

Cr  ■  siurface  tension  of  liquid 


We  define 


-=  4TrfV  (1) 

Then  •=  (Z) 

It  is  seen  that  has  a  maximum  value  at 

9  =  g*  •  (% 

idiere  g*  is  defined  as  a  critical  sized  embryo,  being  the  size  that  involves 
the  greatest  change  in  free  energy,  and  hence,  most  difficult  to  form. 

For  a  stable  vapor  for  which  ( ^  opposed  to  a  supersaturated 
one  ),  the  steady  distribution  of  embryos  is  given  by  the  Boltzman 


distribution: 


We  then  assxime,  that  even  for  the  case  this  same  distribution  may 


be  used,  and  thus  obtain 


-  Ce-^^  fnax 


where  C  roughly  equals  the  total  number  of  molecules  in  the  system. 
From  (1),  (2)  and  (3) 


max 


41  d 
3 


where  r  *  r* 


when  g  »  g* 


5. 


The  di8tri'i>ution  (5)  shovs  Mg*  to  correspond  to  the  niniaos  ralue  of 
Mg.  Use  of  tius  distribution  assunes  that  there  are  no  autual  interactions 
between  the  various  eabzyos,  which  can  therefore  be  treated  as  if  th^  were 
isolated  from  (me  another. 

After  writing  the  appropriate  expressions  for  •  (b  t  in  (2),  differ- 

entiation  at  fixed  (T  »  T^),  with&v4^sO*  assuming  gives 

"  '  6  A 

In  •#  =  2o-0a  ^  2.C-  (6) 

^  r-kT  e^RTjr* 

idiere  Og  *  volume/aolecule  in  liquid  phase 
»  gas  constant 

*  tenqmrature  of  liquid  drop 

From  here  on,  we  further  assume  that: 

1)  All  drops  of  a  certain  size  G  >’  g*  are  eliminated  from  the 

system  and  replaced  by  an  equivalent  nuaber  (GN^)  of  single  molecules. 

Thus  the  number  of  drops  of  any  size  remadn  constant.  Hence,  if 

»  number  of  molecules  eviq>orated  ^om  the  surface  of  a  g-sized 
drop/unit  tine/unit  area 


P  - 


number  of  molecules  striking  unit  drop  surface/unit  time. 


This  assumption  requires  that: 


^  is  given  from  kinetic  theory  as: 


f^(2irmkT) 


and  S  s  surface  area  of  a  g-sized  drop.  For  1>  it  is  a  good  assmp-* 

tion  that  .  ai  S^. 

g*“l  g 

Hence,  substitution  for  N  ,  N  .  from  (5)  in  (7)  gives  the  result: 

O  ft  ^ 

=r  fi^(^T)(<t'e'<t*A+2(rOB/rgj  (9) 

We  see  that,  for  g  =  g* 

0<g<r  P 


6. 


We  now  essuae  that  as  given  in  (9)  holds  even  for  a  non-equilih- 
riia  distribution  of  enibxyo  sizes,  and  proceed  to  determine  the  rate  of 
change  froa  a  non>^quilbxlua  distribution  characterized  by  f  •  Tsking 
account  of  condensation  on,  and  evaporation  from  the  dr(9,  ve  define  as: 

I3  = 

equals  rate  of  fpmation  of  g-sized  drops  due  to  condensation  on  (g-l)  sized 
drops  and  evaporation  of  g-sized  drops.  Siailarly, 


equals  rate  of  decay  of  g-sized  drops  due  to  condensation  on  g-sized  drops 
and  evaporation  of  (g^l)  sized  drops. 

Hence,  net  rate  of  formation  of  g^sized  drops  is  given  from  (10)  and 
(11)  as:  « 

~  V'  (12) 

In  steady  state  ^g/  dt  *  0 . 

For  g^lO,  f  can  be  treated  as  a  function  of  a  continuous  variable, 

o 

g,  and  (Ig  -  replaced  by  -  ii/j*. 

The  appropriate  boundary  condition  is  that  at  g  ■  G,  f(g)  *  0,  and  an 
approximate  solution  of  (12)  can  be  expressed  as 


exp(-4irffr*y3kT) 


where  the  further  assumption  has  been  made  that  I  is  independent  of  g  for 
8<0*»  r*  is  given  by  (6). 

Certain  modifications  to  this  theory  have  recently  been  suggested  by 
(28) 

Lothe  and  Pound'  ^vho  have  pointed  out  additional  terms  which  should  be 
Inlcuded  in  the  free  energy  of  fonaation  of  critical  nuclei  as  mentioned 
above: 


a)  They  believe  that  errors  due  to  assigning  macroscopic  proper¬ 
ties  to  the  embryo  may  be  reduced  by  taking  account  of  the  entropy 
reduction  accompeuiying  the  separation  of  say,  g  molecules  from  a  larger 
liquid  bulk  of  gN  molecules.  This  would  add  to  the  free  energy  of 


fomation  of  critical  sized  eaibzyos,  a  ten  kT  In 

and  lead  to  a  consequent  reduction  in  the  equilihriiam  distribution  of 
critical  nuclei  of  an  order  of  magnitude 

b)  Since  6  degrees  of  freedom  are  required  to  energize  the  eid>ryos, 
conservation  of  degrees  of  freedom  requires  that  6  degrees  of  freedom 
in  the  individual  molecules  be  de-activated.  nils  leads  to  a  further 
Increase  ^  ^  given  by 

“  i*  "Ts 

and  hence,  to  a  decrease  in  the  equilibrium  concentration  of  critical 

-2 

nuclei  by  a  factor  of  about  10  . 

c)  The  embryos  also  have  3  translational  degrees  of  freedom  and  a 

certain  translational  energy  associated  with  these,  which  lead  to  a 
decrease  in  given  by 

=  -  kTlo[(2rmkTf'!a 

where  m  »  molecular,  mass  of  embryo, 

n  *  molecular  volume  of  embryo  in  gaseous  state  at  T. 

For  water  vcpor,  this  is  of  the  order  of  -  2k  kT  when  g*:^100  molecules, 
and  leads  to  an  increase  in  I  by  a  factor  of  about  10^^. 


d)  The  rotational  degrees  of  freedom  of  the  embzyos  also  contrib¬ 
ute  to  the  energy  and  entropy  of  formation  of  critical  nuclei.  Assum¬ 
ing  the  molecules  in  the  embryo  to  be  rigid,  and  the  droplet  spherical, 
so  that  s  Ig  s.  ^  *  moment  of  inertia),  we  obtain  a  further 
contribution  to  given  by 

which  is  roughly  «  -  21  KT  for  water  vapor  droplets  with  g*  *  100  at 

Q 

300  ®K;  and  increases  I  by  a  factor  of  about  10  . 

Sxanming  the  factors  (a)  to  (d),  we  have 

Afj  +  ^  -45kT 


for  voter  droplets,  corresponding  to  an  increase  in  the  distribution  of 
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critical  nticlei  by  o  factor  of  roughly  10  . 


More  recently,  Oriani  and  Sundqiiist^  .  by  similar  considerations  to 

l8 

the  above,  have  arrived  at  a  correction  factor  of  about  10  for  water  vapor; 
^  r>  11  iri -j-  ■«  m  ar  (19).  (29) 

1  M  M  t  A  4  ^  n  ^  n 


and  Courtney 


following  the  'conservation-of«^egree8-of>freedom- 


k  6  9 

approach',  has  come  iq>  with  correction  factors  of  about  10  ,  10  ,  and  10 
for  aonatoaic,  linear  and  non-linear  polyatoeiic  condensing  molecules,  res¬ 


pectively. 


Further  Growth  of  the  Drops 


(5,  10) 


Ibe  growth  rate  of  the  drops  is  obtained  by  considering  mass  transfer 
between  vapor  and  droplets .  Let  ^  ^  the  fraction  of  molecules  striking  a 
dr(9  siurface  that  condense  on  it.  The  mass  flux  from  the  dr(^  equals  meiss 
flux  to  it  from  an  environmental  condition  T^,  P^,  corresponding  to  drop 
teaqierature  and  pressure.  With  defined  by  (8)  the  mass  fltix  from  the 
drop  surface. 


S*»//2TrRT 


where  »  saturation  pressure  corresponding  to  drop  temperatxure ,  T^.  But 
mass  flux  to  the  surface  of  the  drops 

=  .  -j 

/.  B,  continuity:  =  4^rr^|!fc  ~ 


it  “ 


The  drop  teoqperature  is  calculated  by  a  balance  of  energy  gluxes .  We 
assume  an  accommodation  coefficient  of  1,  which  inq>lies  that  all  of  the 
energy  associated  with  molecules  striking  a  drop  surface  is  reflected,  at 
the  tesqperature  corresponding  to  T^.  These  considerations  lead  to  the 


expression 


/  41 


^  /tt 
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It  is  foond  that  dTj^/dt  ^  dr/dt ,  and  hence  ve  assuoe 
Slininating  dr/dt  from  (13)*  (1^)*  plus  use  of  (6) 


II.  EXPERIMEHTAL  WORK 


Tjj»Tjj(r). 

(15) 


A.  Preliminary  Considerations 

It  vas  desirable  for  experimental  piirposes  to  choose  a  test  fluid  which 
condenses  with  marked  property  changes  and  which  would  therefore  render  the 
detection  of  condensation  fairly  easy.  This  purpose  would  be  adequately 
served  by  a  fluid  which  first  has  a  high  nucleation  rate,  and  which,  from 
a  given  stagnation  condition,  approaches  saturation  rapidly  since  it  is 
desirable  to  have  condensation  at  fairly  low  Mach  numbers.  An  e:%amination 
of  the  nucleation  rate  equation  shows  it  to  be  strongly  dependent  on  the 
eiqponential  index,  and  hence,  a  useful  parameter  of  nucleation  rate  might 
be  taken  as  K.  *  •  A  low  value  of  this  parameter  would  indicate  a 

1  ^  R 

very  high  tendency  to  form  condensation  nuclei  very  rapidly.  Ihe  second 
requirement  of  rapid  i[q)proach  to  saturation  is  governed  by  the  parameter 


K, 


Kfo 

m 

■vr 


obtained  from  gas-dynamic  considerations.  A  better  grasp  of 

(17) 


its  actual  significance  is  obtained  from  the  following  illustration 
The  Clausis-Clapeyron  Equation  gives,  along  the  saturation  cuirve: 

(16) 


§■  - 


ar  W 

For  isentropic  expeuision. 


Hence, 


- 

T  /&-I 


= 


^ArT 


(17) 

(18) 


For  high  values  of  /C  T  the  saturation  line  is  much  steeper  than 

I  g  p 

the  isentropGs,  and  hence,  for  decreasing  temperature  along  a  given  isen- 
trope,  rapid  approach  to  the  saturation  line  is  achieved,  lypical  values 
of  the  two  parameters  for  ammonia,  and  a  few  other  vapors  are  listed  below 


for  comparison. 
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Trt)le  I. 


0^ 

c  T 

T  ®8 

L 

P 

Aaannia 

.3 

2.5 

U56 

Freon  12 

1.3 

l.h 

U70 

V 

1 

2.8 

706 

Merciuy 

196 

U.5 

920 

normalized  vlth  H^O  *1 

CoB^arison  of  Condensation  Paraiceters  for  a  Few  Flixids 

Initial  estimates  of  the  boundaxy  layer  thickness  at  the  throat  of  a 
converging  passage,  assuming  a  linearly  increasing  velocity,  and  a  reason¬ 
ably  arbitrary  initial  momentum  thickness,  S  ,  were  made  as  a  guide  to- 
wards  the  selection  of  a  nozzle  geometry  which  would  truly  Justify  neglec¬ 
ting  boundaxy  layer  corrections,  ^e  calculations  show  that  the  asstuned 

6L  is  relatively  insignificant.  The  method  used  was  that  of  Launder^ 

o 

in  which  effects  of  coB^>ressibility  of  the  fluid  were  neglected.  The  cal¬ 
culated  ^*/D  at  the  throat  of  the  nozzle  was  of  the  order  of  1%, 
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B;  Apparatus  and  Experimental  Procedure 

The  cofl^lete  echenatic  of  the  Apparatiis  is  shewn  in  Figure  1,  with  more 
details  shown  in  Figtares  2.  Aanonia  is  taken  from  the  tank.  A,  where  it 
exists  in  equilibriiun  with  its  vapor  at  a  pressure  of  ll^a  psig  and  70^F. 

It  is  then  led  into  the  boiler,  B,  the  flow  being  controlled  by  the  valve 
(I).  The  boiler  consists  of  about  l6  turns  of  5/8“  O.D.  stainless  steel 
tubing  through  which  condensing  steam  is  passed  at  pressures  ranging  from 
10  to  20  psig.  The  ammonia  flows  on  the  outside  of  the  coils.  It  was 
found  that,  with  this  arrangement,  steady  state  was  more  easily  obtained 
than  if  the  flow  paths  of  the  steam  and  ammonia  were  interchanged.  A  sight 
gauge  showed  the  level  of  ammonia  in  the  boiler.  As  much  as  possible,  the 
level  was  kept  the  same  throughout  an  eiqterlmental  run.  The  cyclone-action 
separator,  C,  served  to  remove  all  liquid  droplets  left  in  the  vapor,  after 
leaving  the  boiler.  The  ammonia  then  passes  to  the  test  section  through 
the  needle  valve  (II),  which  was  used  to  control  the  stagnation  pressure. 

The  test  section  consists  of  a  stagnation  tank,  6"  I.D.,  with  three 
perforated  screens  held  in  position  by  sleeves,  and  arranged  as  shown  in 
Figure  2b.  The  tank  is  fed  by  about  l6  turns  of  5/8"  O.D.  aluminum  tubing, 
and  the  entire  assembly  is  immersed  in  a  rectangular  ta^  filled  with  water 
at  various  temperatures.  The  stagnation  tank  carries  a  pressure  gauge,  a 
thermometer  and  thermocouple.  From  the  settling  tank,  the  ammonia  goes  into 
the  test  nozzle,  along  which  measurements  of  pressure  distributions  are  made. 

Two  nozzles  have  been  used  in  the  tests,  an  axi-symmetric  nozzle,  shown 
in  Figure  2b,  and  later  a  two-dimensional  one.  The  one -dimensional  design 
was  based  on  the  voi-..  of  Wang  and  Co,  who  described  a  method  of  obtaining 
uniform  flows  at  the  outlet  of  a  converging  chemnei;  with  a  view  towards 
having  a  uniform  flow  at  the  throat.  The  diverging  part  of  the  nozzle  was 
conical,  although  this  meant  having  a  discontinuity  of  curvature  at  the 
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throat.  It  vas  felt  hoire'ver  that  this  could  not  introduce  any  great  disturb'- 
ances  to  the  aeasured  pressure  distributions.  Pressure  taps,  20/1000”  I.D. 
vere  located  at- 1/1(''  apart  along  the  vails  in  the  diverging  section .  The 
diverging  part  vas  3”  long,  had  a  throat  diameter  of  2/10”  and  an  exit 
diameter  of  27/100” «  Such  small  dimensions  vere  necessitated  by  mass  flov 
considerations.  Itofortxmately,  hovever,  this  nozzle  did  not  prove  very 
useful.  First,  for  a  nozzle  of  such  small  dimensions,  the  pressure  taps 
have  to  be  perfectly  smooth  and  rounded  off.  ^Any  small  burrs  or  irreg¬ 
ularities  in  the  taps  proved  to  have  a  considerable  influence  on  the  pressure 
distributions.  Hence,  after  any  one  run,  the  deposited  moisture  usually  had 
a  disturbing  effect  on  the  taps,  and  the  follovlng  set  of  pressure  measure¬ 
ments  usually  differed  considerably  from  the  previous  one,  enplpying  al/most 
Identical  tq>stream  conditions .  An  attetQ>t  to  solve  this  problem  vas  made 
by  reaming  out  the  nozzle  after  every  fev  runs>  but  this  in  itself  vas 
self-defeating,  since  it  meant  that  the  runs  vere  necessarily  irreproducible 
becaiise  of  the  constantly  changing  geometry.  It  vas  to  resolve  this  dilemma 
that  the  tvo-dlmenslonal  nozzle  vas  constructed.  This  is  shovn  in  Figure  3. 

The  inlet  consists  of  circular  arcs,  vith  a  1”  radius.  The  angle  of 
divergence  can  be  varied  betveen  about  1^  to  U**  total  included  angle* 

Pressure  taps,  16/1000”  diameter  vere  spaced  1/U”  apart  near  the  throat, 
and  1/2”  apart  further  dovnstream.  The  throat  vidth  vas  about  15/100” 
independent  of  divergence  angle,  and  the  nozzle  height,  uniform  everyvhere, 
vas  1/5".  The  nozzle  can  be  taken  apart  and  cleaned  betveen  runs,  a  great 
advantage  over  the  one-dimensional  one. 

The  annonia  is  finally  discharged  by  dissolving  it  in  a  steady  stream 
of  cold  vater  and  leading  the  mixture  into  the  drain.  The  air  passages 
(Ref.  Fig  la)  are  provided  for  bloving  out  either  the  test  section  or  the 
entire  apparatus. 
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!nie  procedxure  for  obtaining  a  con^lete  set  of  preeanre  distribixtions 
is  as  foUovs:  (ref.  Fig.  la).  The  water  in  the  rectangular  tank  is  kept 
at  the  appropriate  teBQ>erature .  Steam  is  turned  on  to  a  suitable  pressure 
level,  usuadly  about  10  psig,  and  left  tunning  for  about  5  minutes  to  warn 
up  the  coils  in  the  boiler.  The  cold  water  is  next  turned  on  and  a  steady 
stream  maintained,  with  a  little  overflow  from  the  qiienching  tank,  D.  Next, 
the  outlet  valve  from  the  aanonia  tank  is  opened,  so'  that  valve  I  experi¬ 
ences  the  full  tank,  pressure .  With  the  needle  valve  II  now  slightly  opened, 
valve  I  is  opened  fairly  wide,  and  immediately,  the  exhaust  valve  III  is 
opened.  Fine  regulation  of  the  needle  valve  now  gives  any  desired  stag¬ 
nation  pressure,  a  rough  value  of  which  is  read  on  the  pressure  gauge.  The 
valve  III  is  necessazy  to  prevent  suck  back  of  water  from  the  quenching 
tank  when  the  system  is  shut  off.  The  stagnation  teaq>erature  can  be  varied 
between  two  extremes,  neamely,  the  triple  point  of- water  and  the  boiling 
point.  Initially,  the  water  in  the  tank  is  heated  electrically  to  near 
boiling,  and  for  subsequent  runs,  varying  intermediate  temperatures  are 
obtrined  by  emptying  differing  amoxmts  of  the  hot  water,  and  refilling  the 
tank  with  cold  water.  Eventually,  the  water  is  cooled  to  normal  cold  water 
temperature,  from  which  point  its  temperature  can  be  lowered  still  further 
by  adding  varying  amounts  of  ice.  Strictly  speaking,  the  intermediate 
temperatures  cannot  be  held  constant,  but  since  it  takes  less  than  2  minutes 
to  take  a  set  of  pressure  distributions,  it  is  assumed  that  the  tank  temper¬ 
ature  cannot  vary  much  in  this  interval.  This  assumption  is  borne  out  by 
the  experimental  data. 

The  pressure  readings  are  taken  with  a  transducer,  utilising  a  rotary 
switch  device  shown  schematically  in  Figure  4.  The  switch  consists  of  two 
discs,  the  top  one  rotating  on  nylon  bearings  mounted  axially  through  the 
bottom  one.  The  latter  has  2k  holes  drilled  through  it  circumferentially. 
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and  each  connected  to  a  pressure  ti^  by  flexible  tubing.  Each  pressure  hole 
in  the  lover  disc  is  individually  sealed  by  an  0>ring  between  the  two  discs, 
and  a  single  search  tube  through  the  t^per  disc  picks  up  the  pressure  in 
each  tap  as  it  turns  around,  and  transmits  it  to  the  transducer.  Callibra- 
tioo  shows  the  transducer  response  to  be  linear  with  pressxire.  The  signal 
is  fed  into  an  aaqplifier  with  a  suitable  outlet  for  connection  to  an  x-y 
recorder.  The  recorder  is  calibrated  to  record  the  numerical  valiies  of 
the  pressures  directly. 

The  stagnation  teng>erature  is  read  simultaneously  with  a  thermometer 
and  a  Chromel-Alumel  thermocouple,  both  of  which  have  been  previously  cali¬ 
brated  and  found  to  check  with  each  other  to  within  1/2^F. 

The  average  mass  flow  in  most  of  the  runs  was  about  3  cu  ft/hr,  with 
stagnation  pressures  varying  between  30  and  70  psia,  and  stagnation  temper¬ 
atures  between  32°  and  190°F. 

C.  Precision  of  the  Measurements 

The  calibration  of  the  transducer  showed  it  to  be- linear  with  pressure, 
as  illustrated  in  Figiure  5*  To  allow  for  the  possibility  that  the  straight 
line  might  shift  parallel  to  Itself,  the  zero  pressure  point  was  usuedly 
checked  before  each  set  of  readings.  It  is  assumed  that  the  gedn  of  the 
d.c.  aBg>lifier  does  not  depend  on  the  magnitude  of  the  excitation  voltage, 
and  hence,  a  linear  input  would  also  give  a  linear  output.  However,  since 
the  calibration  of  the  amplifier  output  was  achieved  by  setting  the  dial 
pointer  to  correspond  to  a  known  pressure  from  a  mercury  manometer,  there 
exists  two  sources  of  uncertainty;  first,  in  reading  the  nanometer  levels, 
and  secondly,  in  setting  the  pointer  to  correspond  to  the  known  pressure. 
These  two  errors  could  of  course  be  in  opposite  directions,  and  hence  can¬ 
cel  out,  but  assuming  they  are  in  the  sane  direction,  and  allowing  an  \incer- 
tainty  of  t  O.l*"  for  the  mercury  levels,  the  cumulative  error  could  be  of 
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the  order  of  i  0.7  psia.  There  Is  a  third  error  InTolred  in  reading  the 
recorded  pressures,  which  can  only  be  read  conrenientlj  to  about  ±0.2  psi. 
Thus,  the  worst  possible  uncertainty  in  the  pressures  shown  on  the  data  aay 
be  assuned  to  be  not  greater  than  0.9  or  about  2%»  ■  Concerning  the  speed 
of  response  of  the  entire  sensing-aa^lificatlon-recording  systen,  suffic> 
lent  tiB«  was  generally  allowed  for  a  flat  pressure  response  to  be  obtained 
at  a  given  pressure  tap  before  switching  on  to  the  next  pressure  tap.  The 
typical  dturation  of  tine  spent  at  a  presstare  location  was  sbout  3  seconds; 
during  which  the  pressure  was  observed  to  attain  a  flat  steady  value  within 
the  first  3/k  of  a  second. 

The  accuracies  regarding  the  teniierature  neasureiaents  are  perhaps  a 
little  nore  precise.  It  is  difficult  to  balance  the  potent! osieter  any 
closer  than  to  within  t  0.02mv,  or  about  t  l^F.  An  additional  error  of 
nay  be  assuaed  due  to  the  reference  Jxmction  not  being  exactly 
at  32*’F,  this  being  alwsys  in  the  positive  direction.  Except  in  a  few  runs, 
tenperatxures  at  the  beginning  and  at  the  end  of  a  run  did  not  differ  by 
nore  than  1.5*’F  and  most  of  this  is  estimated  to  result  from  the  error  of 
balancing  as  mentioned  earlier.  The  temperatures  as  given  in  the  data  may 
therefore  be  assuaed  to  be  accurate  to  within  -1.0  to  ♦  1.5®F,  A  further 
point  may  be  made  as  regards  the  purity  of  the  ammonia.  The  commercial 
grade  that  was  used  was  certified  by  the  manufacturer  to  be  of  99*95%  pur¬ 
ity,  contcdning  less  than  50  ppm  of  water  vapor,  and  less  than  5  Pl«  of  oils 
and  other  non-condensibles ^ .  The  nucleation  that  led  to  the  condensation 
may  therefore  be  assumed  homogeneous  as  explained  earlier. 

D.  Experimental  Results  and  Discussion 

The  upper  curves  in  Figures  6  and  7  show  the  experimental  pressure 
distributions  obtained  with  the  axisymmetric  nozzle.  The  lover  curves  are, 
respectively,  the  theoretical  pressure  distribution,  and  the  isentropes. 
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Figures  8-13  shov  the  pressure  distributioos  from  the  tvo-dimensional 
nozzle.  The  curves  in  each  figure  correspond  approximately  to  the  same 
stagnation  pressure,  but  vith  differing  stagnation  teiq>erature8 .  They  were 
all  obtained  vith  a  total  nozzle  included  angle  of  1^  •  The  lowest  curve 
corresponds  to  a  high  enough  8\q>erheat  of  the  vapor  so  that  no  condensation 
occurs  in  the  nozzle.  Ve  refer  to  this  as  the  ‘nou-condensing*  curve,  and 
use  it  as  a  reference  for  the  other  curves.  Queaitatively,  Figures  8-15 
are  vexy  similar,  and  hence,  the  following  remark  .  referred  specifically 
to  Figure  10,  is  applicable  to  all  of  them 

1)  The  non-condensing  curve  lies  in  general  above  the  isentropic 
curve  (not  shown)  corresponding  to  the  geometry  of  the  nozzle.  This 
is  mostly  due  to  frictional  effects,  resulting  in  an  effective  nozzle 
area  different  tron  the  geometric  area  at  any  cross-section.  The 
fX'lctional  effects  are  more  pronounced  in  the  supersc  lic  section  of 
the  nozzle  than  in  the  subsonic  part.  For  the  same  reason,  the  meas¬ 
ured  pressure  ratio  at  the  geometric  throat  is  much  greater  than  the 
e^cted  isentropic  value  of  0.5^^»  for  an  Isentropic  exponent  of 

k  «  1.3,  shoving  that  the  effective  throat  is  actually  downstream  of 
the  geometric  throat. 

2)  As  the  stagnation  temperature  is  sufficiently  lowered  below 
the  value  corresponding  to  the  lowest  curve,  the  measured  p/p^  at  some 
point  in  the  nozzle  becomes  higher  than  the  corresponding  p/p^  in  the 
reference  curve.  Downstream  of  this  point,  measured  p/p^  are  every¬ 
where  higher  than  in  the  reference  curve.  Ibis  deviation  of  measured 
pressure  ratios  from  the  reference  values  is  taken  as  a  sign  of  con¬ 
densation  from  the  vapor.  The  condensed  fluid  releases  its  enthalpy 
of  vaporization  to  the  remaining  gaseous  phase,  which  results  in  a 
change  in  stagnation  temperature  of  the  fluid.  From  one-dimensional 
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considerations  and  for  effect  of  heat  addition  alone,  it  is  known  that 

Jp-  _  _JcMf  ,  , 

Hence,  in  si^erscmic  flow,  such  a  heat  release  to  the  wapor  asy  he 
expected  to  lead  to  a  pressure  rise  as  shown  hj  the  data. 

3)  As  the  te^peratiure  is  further  lowered,  the' point  of  pressure 
dewiation  aowes  iq>8treaai,  and  the  deviations  firoa  tke  non-condensing 
curve  at  ansr  cross-section  becones  greater,  as  seen  ttm.  equation  (19), 
when  M  ->1  frosi  values^  1. 

k)  With  further  decrease  in  stagnation  tea^erature  the  pressure 
curve  actually  begins  to  exhibit  a  definite  adnianm  and  aaxlaun,  and 
the  location  of  this  alninum  is  in  general  different  froa-  the  point  of 
deviation  of  the  curve  from  the  non-condensing  ciurve.  We  have  used 

the  convention  that  if  the  distribution  does  not  show  a  definite  lain- 

*  ' 

iamn  and  amxlaniBgi,  then  ,the  mset  of  condensation  is  defined  as  the  point 
of  departure  fron  the  non-condensing  curve.  Otherwise,  the  adninna 
point  of  the  curve  is  chosen. 

It  has  been  observed  that  the  non-condensing  curves  sure  thiwselves 
dependent  on  the  pressure  level.  This  effect  is  illustrated  in  Figure  l6. 

At  the  lowest'  stagnation  pressiure  the  non-condensing  curve  follows  the  path 
A  -  C.  With  increasing  stagnation  pressures,  however,  it  begins  to  approach 
the  path  B  downstre.’jn  of  the  nozzle,  while  still  keeping  vezy  close  to  path 
A  upstream,  until,  at  about  60  psi  stagnation  pressure,  it  actually  does 
follow  the  path  A  -  B.  This  Justifies  the  statement  made  earlier  that 
frictional  effects  are  mo'^tly  confined  to  the  section  of  the  nozzle  down¬ 
stream  of  the  throat.  It  was  further  observed,  however,  that  with  pressures 
between  60  -  70  psi  (Figures  ll»  &  15)  the  non-condensing  curve  shows  a 

t  . 

shift  to  path  D  in  the  upstream  side  of  the  nozzle,  and  follows  path  D  -  B. 

It  is  believed  that  this  further  shift  is  not  due  to  friction  for,  although 
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ve  have  seen  that  ftrictional  effects  are  much  greater  in- the  stq;>ersonic 
than  in  the  subsonic  region,  there  is  hardly  any  noticeable  shift  in  the 
si^rsonic  section  within  the  same  interval  of  pressure-  increase  as  the 
subsonic  path  shifts  from  A  to  D.  Furtheraore,  the  flow -along  path  D  was 
observed  to  be  very  unsteady,  particxtlarly  within  a  saall  distance  of  the 
geosKtric  throat.  This  nay  be  seen  troa  Figure  l8.  As  this  observation 
seens  to  be  independent  of  the  initial  s!;9erheat  of  the-  vapor,  it  cannot 
possibly  be  due  to  condensation  occurring  too  near  the  throat,  as  we  had 
thought  earlier.  The  xmsteadlness  is  in  fact  greater  for  the  higher  super¬ 
heats.  It  is  postulated  that  these  disttirbances  could  have  been  caused  ly 
non-unifomities  in  the  entering  flow  which,,  iwhile  unhotlceable  at  lower  pres¬ 
sures,  became  magnified  as  the  pressure  increased. 

Figure  17  shows  pressure  distributions  obtained  frcn  the  two-dimen¬ 
sional  noxzle  with  a  total  included  angle  of  3°.  They  illustrate  the  same 
qualitative  feattures  as  those  earlier  described.  The  initial  superheats 

and  degree  of  supersaturation  achie-ved  Just  before  condensation  are  shown 

(21) 

on  the  charts.  Contrary  to  the  findings  of  Duff'  'it  does  not  seem  that 

there  is  any  definite  correlation  between  the  initial  superheat  and  final 

degree  of  si;q>ersaturation  attained,  although  at  the  higher  pressures  used 

in  our  experiments,  the  supersaturation  tended  to  increase  as  the  shock 

moves  towards  the  throat,  but  even  this  is  contrary  to  Duff's  results.  In 

this  respect,  our  present  findings  cuw  in  egreement  with  the  experiments  of 
(9) 

Blnnie  and  Woods  on  steam. 

It  may  be  asked  at  this  point  what  effect  the  wall  has  on  the  conden¬ 
sation?  A  rough  heat  transfer  analysis  showed  that  the  nozzle  wall  is 
superheated  relative  to  the  free  stream  pressure  and,  hence,  cannot  pos¬ 
sibly  initiate  the  condensation  process,  which  may  thus  be  regarded  as 
essentially  a  free-stream  phenomenon. 
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III.  NIMERICAL  CALCULATIONS 
A.  Isentroplc  Expansion 

Before  condensation  occurs*  the  flow  through  the  nozzle  is  considered 
to  be: 

one-diiaensional,  steady*  Isentroplc. 

The  stagnation  conditions  in  the  stagnation  tank  before  the  entrance  of  the 
nozzle  are: 


P  ■  Pq,  T 


T^*  V 


V  ■  0. 

O 


The  cross-sectional  area  at  any  place  x*  A  ■  f(x)t  is  given  by  the  geometric 
design  of  the  nozzle.  From  this  A  ■  f(x),  the  unknown  properties  p,  T*  V 
■  f(A)  or  f(x)  at  any  point  x  along  the  nozzle,  can  be  calculated  from  the 


equations  for  choking  flov,  (see  Ref.  (22): 


Jl  ^  -L  r  ^  llic^ 

A*  M  L  K+i  J  K 

[I- 


-B.  * 
Po 


[ 


y  (20) 


i 


M 


With  the  Mach  number  M  herefrcii.  '-^.nown,  the  velocity  V  is  given  by: 

V  ■  M  .^k  g  R  -  M  •  a,  where  a  ■  y  k  g  R  T  ^ 
is  the  local  velocity  of  sound  in  the  vapor. 

The  mass  flow  for  choking  flow  at  the  throat  is  given  by: 

where  lu  ^ 

/iHWl  (  K.'h  f  J  Ki"! 

where  k  ■  c^/c^  is  the  isentroplc  exponent  at  the  conditions  at  the  throat. 
B.  Nucleation  Rate  I 

The  expansion  is  considered  to  remain  isentroplc  until  the  nucleation 

rate  I  ■  f(p*  T)  reaches  a  characteristic  critical  value  I  which  we  chose 

cr 

to  be  10^®  nuclei^t?  sec).  The  onset  of  condensation  was  found  to  be  rel¬ 
atively  independent  of  the  assumed  value  of  I  when  I  wou  varied  netween 

cr  cr 

the  limits  from  1  to  10^°  nuclei^^  sec}.  As  soon  as  I  *  f(p,  T)  ^  I  ,  the 


cr' 
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isentropic  cfld.culati<»i  is  finished  and  condensation  effects  are  taken  into 
the  further  calculation  procedure. 


The  egression  for  the  nucleation  rate  I 


3kT 


vhere  ^  «  r — ^  77. .  ■  .  '  ' 

^L'RT-ei^rypg^ 

is  the  critical- dr<q;>  radius  fron  equation  (6).  Drops  greater  than  r*  can 

grow,  drops  smaller  than  r*  vaporize  again.  The  surface  tension  6*  is,  since 

the  droplets  ere  very  small,  assumed  to  be  considerably  different  froai  the 

flat  film  surface  tension  6*^  .  A  brief  discussion  about  necessary  cor> 

rections  to  the  surface  tension  follows  in  Part  IV,  B.  Here  the  assximption 

(26 

according  to  Tolman  is  tised: 

where  is  a  length  between  0.50  and  1*20  of  the  molecular  radius  r^. 

The  .urfwe  teneion  of  the  flat  fll«  eurfwie  ^  is  dependent  on  the 
teoperature  T,  see  Ref*  (2U): 

/-  Ke  CtJ  -T) 


S'*,* 


^  'If 

■” . . 

*>  Eotvos  constant 

T^^Tj^  *»  critical  temperature 


where 


A  s  moleculcu:  weight 
*  liquid  density 

Nov  the  quantity  y,  is  defined  so  that  r*  «  «r/y  ,  where 


then  the  surface  tension  becomes 

6*  •  6^  —  *= 

Substitution  of  ^  and  y  into  equation  (21)  gives:  i 


t^itution  of  ^  and  y  into  equation  (21)  gives:  rjfeZL  6*'»7 

r  -e'®" 
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The  Molecular  oass  ■  and  the  nolecular  voluac  y  are: 

idiere  the  aolecule  is  assumed  to  be  a  sphere  of  radius  r^; 

-  o,Q035I-^ 

Figure  19  shews  lines  of  constant  nucleation  rate  in  a  p,  t<>diagraB, 
obtained  from  equation  (23) •  It  illustrates  hov  rapidiy  the  nucleation 
rate  increases  during  the  expansion  after  saturation  is  reached.  A  fev  is- 
entropic  e:q>ansion  lines  are  drawn  to  show  this. 

C.  Calculation  of  Pressure,  Temperature  and  Velocity  along  the  Woxzle 
under  Consideration  of  the :Presence  of  Liquid  Parts 

The  following  further  as8uiiq>tion8  are  made: 

a)  one-dimensional,  steady  flow 

b)  no  heat  transfer  frcmi  the  surrounding  wadis 

c)  no  friction  at  the  walls 

For  the  flow  in  the  nozzle  the  basic  equations  are  written  for  a  con¬ 
trol  volume  of  the  length  dx: 


t  Jv  ^ 
Y  ^ 


(26a) 


(26b) 


(27a) 


(27b) 


a)  Continuity  Equation 

am*  const  (26a) 

4- +  X  £^t  +  J! _ (261,) 

Y  Q  e/x  fix  A  cR 

b)  Momentum  Equation 

—A  -  A  (27a) 

X  '^1'  +  '  .X. 4^  (276) 

c)  Energy  Equation 

+\olVmO  (28«) 

(I  -f-)  -vf -0 

d)  Equation  of  State  of  a  Perfect  Gas 

p  m  ^  R  T  (29) 


(29) 
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Sy  definition  i 


V 

iicaRT' 


and  K 


Si 

Cy 


We  introduce  a  diisenaionleas  parameter 

^ 

After  zvplacing  ^  by  p  and  T  and  introducing  the  Mach  number  M  the  equations 
(26),  (27)  and  (28)  becoone: 


Continuity:  /  c/V  ^J_c/T  _  ^J_dA  __/ _ 

V  efx  p  tU  tW  ~  A  W  (/-.(A)  oix 


Momentum: 


Energy: 


>(30) 


Remark:  An  intermediate  step  in  the  energy  equation  gives: 

eU  ‘  t  0/x  j  dx  ^  dx 


ctx 

vhere  the  second  tern  is  neglected  because  of  the  smedl  amount  of 


From  these  the  three  unknown  properties  p»  T,  V  can  be  estimated: 

tlx  0-t*)^  A^j 

M.  s  ^  SL  ..Lels.  /(3i) 

dx  Kh^  j. 

■Sx  “  /-JT-  -p- w  +  ^  -sr*/ 

In  order  to  find  the  distribution  p,  V,  T  »  f(x)  the  derivative  (^dx  has 
to  be  determined.  (A  and  dA/dx  are  known  from  the  geometry  of  the  nozzle). 


D.  The  Liquid  Part  M  and  Its  Derivatives 

The  liquid-mixture  mass  ratio^  is  given  at  any  instantaneous  point  x 
along  the  nozzle  by  the  lolloving  considerations: 

1)  small  droplets  (nuclei)  of  radius  r^  have  been  built  on  a 
nucleation  rate  I  *  f(p,  T)  in  the  past  length  along  the  nozzle. 

(  O-^J  -X). 

2)  those  nuclei  grew,  and  their  radius  r  increased  by^r/dx  » 

f(p,  T)  steadily  from  the  point  where  they  were  built  to  the 
point  ->  x) . 
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It  is  assumed,  that  at  any  point  x  all  the  dr<9lets  grow  at  the  same  rate: 

X  dr/dx,  regardless  of  their  size. 

The  experessicn  for^  at  the  point  x  is  given  by  the  following  expres¬ 
sion,  where  ^  and  ^  are  dummy  variables,  as  shown  in  the  sketch: 


After  this  ve  have  the  following  system  of  necessary  equations: 


2k, 


drt 

* 

elB 

Hx 

id 

ttx 

six 

iB- 

dX 

dV 

W 

etr 

dx 

E.  ' 

-h  IAIk)  •  Bn  -  W 
2  +JA(x)  •■J'iA) 

O.  /c/f*  r/. _ /_  _±f^i7 

^  -I  L(  ^  A  ttxj 

_ )L.  .J_et 

kM*-  P  qIx 


/r-/ 


According  to  Ref.  (13)  the  radius  of  the  newly  built  droplets  r^  is 


about  1.3  times  the  critical  radius  r*: 

^  ^  i  «» 

%^/.3  T  -  -ft-Tfrrzc 


The  growing  rate  of  the  drop  radius  dr/dx  is  obtained  from  the  equations 


(13)  and  {Ik): 

c/r  _  3i 


[rt,,-rt] 


where  the  drop  temperature  Tj^  is  given  by  equation  (15): 

In  this  expression  the  drop  radius  is  xinknown.  Here  it  is  assumed  that 
is  an  average  drop  radius  of  all  the  drops  which  strike  the  cross-section 
area  A  at  the  point  x.  The  total  surface  S  of  all  these  drops  passing  the 
area  A  at  the  point  x  per  second  is  gi  ren  by  equation  (33b): 


the  point  x  per  second  is  gi  ren  bj 

3  JIA(0  clf 


and  the  total  mass  of  all  those  droplets  is  given  by  equation  (32): 

So  the  avereige  drop  radixis  becomes 

rzx  ^  3 


u. 


t 

ion 
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The  drop  tentperatiure  has  to  be  conputed  from  equation  (15).  Since 
drop  *'**‘®^®  saturation  pressure  corresponding  to  Tjj,  and  all 

the  properties  functions  of  T^,  it  is  cooplicated  to  solve 

analytically  for  T^.  (Besides  thic^the  surface  tension  S'  depends  on  the 
average  drop  radius  ^  too.)  Iberefore,  equation  (15)  is  solved  iteratively 
starting  with  a  first  guess  *  T. 

F.  Hunerical  Solution  of  the  Differential  Equations  (37) 

Equation  (37)  is  a  systun  of  seven  simultaneous  first  order  differen¬ 
tial  equations.  They  are  solved  step-vise  by  a  fourth  order  Runge-Kutta 
method,  which  is  ready  to  use  and  written  up  as  the  Subroutine  #1381  "RURGE" 
of  the  M.  I.  T.  Coaputation  Center.  This  program  adjusts  the  initially  given 
step-size  according  to  the  demanded  permissible  relative  errors,  which  we 

took  as  0.00001.  The  starting  values  are  given  at  the  point  x,  where  the 

vo; 

nucleatlon  rate  I  ^  I  is  reached.  There,  the  starting  values; 

cr  I 

*1(0)  *  *(o)  *  ®(o)  ■  ^*(0)  *  O'  '“(o)-  ''(o)>  ■'(o) 

are  the  properties  which  are  calculated  by  the  isentrppic  expansion. 

Remarks;  The  physical  and  thermodynamic  properties  for  Ammonia  were  taken 

from  the  following  references: 

Surface  tension  :  References  (2l|)  and  (25) 

Saturation  pressure  Pg^ 

Liquid  density  ^ 

Specific  heat  c^ 

Vaporization  heat  h^^ 

IV.  VARIATION  OF  THE  CALCULATIC^  PROGRAM 
A.  Correction  for  Frictional  Effects 

Because  of  the  discrepancy  between  the  non-condensing  curve  and  the 
isentropic  line,  based  on  the  geometry  of  the  nozzle,  am  effective  nozzle 
profile  has  been  selected  based  on  the  meastired  pressure  distribution: 


^  Reference  (23) 


The  relation  between  the  neasxured  p/p^  and  the  effective  A/A*  is  given 


*  w"  f 


K-i 


2 

K  +  l  J 


figure  20  shows  a  coo^arison  between  the  isentropic  pressure  distrib¬ 
ution  for  the  geanetrlcal  nozzle  profile,  and  that  corresponding  to  the 
effective  nozzle  geometry* 

B.  Uncertainties  in  the  Physical  Properties  Entering  the  Nucleation 


and  Drop- Growth  Theo 


e)  Surface  Tension  6*  : 

The  correction  to  the  flat  film  surface  tension,  used  in  the  calcul¬ 
ation  procedure,  was  given  in  equation  (23): 


Herein  the  correction 


to  lie  between  0*50  and  1.20. 


factor  accordi: 


according  to  Tolman^  ^is  recommended 


shews  the  effect  of  the  surface 


tension  on  the  pressure  distribution. 


b)  The  Factor  j  in  the  Drop  Growing  Equations : 

The  factor^  was  defined  ss  the  number  of  condensed  particles  to  the 
total  number  of  arriving  particles  striking  the  drop  surface  area. 

Calculations  were  done  with  constant  values  for  for  different 

values  of  ^  ^  was  modified  to  be  O.OU,  0.20,  0.50  and  0.80.  It 

turned  out  that  the  p/p^-values  for  a  given  x  did  not  change  by  more  than 
1.5  per  cent  in  this  range  ^  »  therefore  the  conclusion  that  the  factor 
quite  insignificant  to  the  results  is  obvious.  For  the  further  cal¬ 


culations 


»  0.50  was  vised  always. 
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c)  A  Coaplete  Change  in  the  Nucleatlon  Rate 

(28) 

Following  the  suggestion  of  Lothe  and  Pound^  ^tbe  nueleation  rate,  I, 

as  calculated  in  the  described  manner,  according  to  equation  (23  )•  was  jsulti- 

17 

plied  by  a  factor  of  10  : 

I 


10” .  I 


eq.(25) 


Thereby  no  correct im  for  the  surface  tension  was  made: 

M 

ISie  results  of  this  change  in  the  nucleation  rate  are  shown  in  Figures 


21b  and  22. 

V.  COMPARISOS  OF  EXPERIMENTS  AND  THEORY 

Figure  20  shows  the  non-condensing  line  and  the  experimentally  measured 
pressure  distribution  for  the  stagnation  conditions  p^  *  30.3  psia;  t^  «  38^F. 
In  the  same  graph  the  theoretical  isentropic  line,  based  on  the  geometric 
design  of  the  nozzle  is  shown.  The  theoretical  pressure  distributions  for 
both  the  effective  areas  from  the  experimental  non-condensing  line  and  for 
the  geometric  areas  are  shown.  One  notices  that  the  pressure  vcdues  for. 
the  two  different  areas  are  considerably  different,  and  since  the  method  of 
using  the  effective  areas  from  the  measured  non-condensing  line  approaches 
the  experimental  results  in  the  order  of  magnitude  much  better,  this  method 
is  further  used  for  all  calculations*  especially  the  difference  in  the  loca¬ 
tion  of  the  throat  is  evident.  The  theoretical  calculations  in  Fig.  20  are 
done  with  if/  ■  0.80  and  ^  *  0.50,  with  no  other  correction  to  the  nuc¬ 
leation  rate,  I,  as  given  in  equation  (23).  It  is  seen  that  the  theory 
shows  a  considerable  delay  in  condensation  as  conqpared  to  the  experiment, 
a  deli^  which  might  have  to  do  with  uncertainties  in  the  surface  tension. 

Figure  21a  shows  how  the  quite  uncertain  correction  term  for  the  sur¬ 
face  tension  /r^  affects  the  theoretical  pressure  distribution.  Based 
on  the  areas  conq>uted  from  the  non-condensing  line  for  ^5  psia,  pres¬ 

sure  distributions  for  the  following  stagnation  conditions  are  shown: 


^o. 

p  •  ^64,  9  psia;  t  *  Uo®F 
o  o 

U7.  0  psia 

1*14.  8  psia  67®F 

For  these  stagnation  conditions  the  theoretical  curves  vith(^  * 
0.30,  O.dO  and  1.20  ajre  drawn. 

The  result  of  this  set  of  curves  shows  that  an  increase  of  <f  /r^,  which 
means  a  decrease  in  the  surface  tension  6^,  moves  the  onset  of  condensation 
towards  the  upstream  direction  of  the  nozzle.  For  the  case  of  9 

psia,  ttfl»0°F  the  value  of  «  1.20  would  predict  condensation  somewhat 

O  M 

upstream  of  the  throat  in  the  subsonic  part  of  the  nozzle,  and  for  this 
case  the  theoretical  calculation  cannot  be  performed  since  the  mass  flow 
confutation  is  based  on  choking  flow  of  the  imcondensed  vepor.  The  experi¬ 
mental  pressure  distributions  are  plotted  into  the  same  figure.  It  ap¬ 
pears  that  the  higher  values  of  -  O.8O  or  1.20  agree  better  with  the 

esqperiments,  showing  that  a  possible  way  to  inf  rove  the  theory  would  con¬ 
sist  in  studying  more  closely  the  influence  of  the  small  size  of  the  drops 
on  their  surface  tension. 

Figure  21b  shows  the  effect  of  multiplying  the  nucleation  rate  by  a 
17 

factor  of  10  .  It  is  seen  that  this  also  brings  the  theoretical  pressure 

distribution  into  better  agreement  with  the  experiment il  distribution. 

In  Figure  22  the  pressure  and  temperature  values  for  the  onset  points 
of  the  condensation  are  shown.  Here,  the  onset  of  condensation  is  sup¬ 
posed  to  be  the  point  in  the  pressure  distribution  curve  along  the  nozzle, 

where  dp/dx  »  0  in  case  of  an  observed  minimum  point,  (see  case  A  in  the 

2  2 

sketch)  or  where  d  p/dx  «  0  in  case  one  cannot  find  a  definite  minimum 
point,  (see  case  B  in  the  sketch). 

The  theoretically  coniputed  pressure  distribution  curves  always  show  a 
definite  minimum  pressure  point,  and  the  locationsof  the  theoretically  pre¬ 
dicted  points  of  condensation  onset  are  drawn  as  continuous  lines  in  the 


29. 


graph.  Three  theoretical  curves  have  been  shown,  two  for  S  /r^  *  0  and 
<$/r„  *  0.80,  with  no  other  correction  to  I,  and  one  in  which  I  has  been 
multiplied  by  10^^,  with  no  correction  for  6* • 


VI.  CONCLUSIOW 

From  the  foregoing^ the  following  conclusions  may  be  made: 

1)  Amionia  shows  less  supersaturation  than  is  predicted  by  the  Becker- 
Doring  (Frenkel-Zeldovich)  Nucleation  Rate  Equation,  the  degree  of  super¬ 
saturation  observed  in  all  the  tests  lying  between  1.5  and  2.5. 

2)  Corrections  to  the  surface  tension  msy  be  a  very  important  part  of 
the  nucleation  theory,  which  therefore  ought  to  be  given  greater  consider¬ 
ation  than  heretofore. 

3)  The  closer  agreement  obtained  between  theory  and  experiment  when 

17 

the  nucleation  rate  equation  is  multiplied  by  10  lends  some  support  to 

(26) 

the  gassiflcation  concept,  proposed  by  Lothe  and  Pound^  ,  and  by  Oriani 
(30 ) 

and  Sundquist  ,  although  it  remains  doubtful  how  much  of  this  correction 
factor  is  due  to  uncertainties  in  surface  tension,  and  how  much  is  due  to 
gassiflcation. 

U)  A  lot  of  scatter  is  evident  on  the  curve  in  which  an  attempt  has 
been  made  to  plot  a  condensation  line.  This  is  mostly  due  to  uncertainties 
in  the  exact  shape  of  the  pressure  discontinuity,  as  may  be  seen  from  the 
curves . 

5)  There  does  not  seem  to  be  any  definite  correlation  between  the 
initial  superheat  and  final  degree  of  supersaturation  attained  before  con¬ 
densation  occuri),  which  may  be  partly  due  to  the  scatter  in  the  results. 


30 


tut  in  the  higher  pressures  used  in  the  eiqaerinents ,  the  8Tq>ersaturation 

tended  to  increase  as  the  pressure  discontinuity  noves  tq^wtresa.  This 

(21) 

finding  is  contrary  to  that  of  Duff  althoxigh  in  agreeaent  with  those  of 
(9) 

Binnie  and  Woods  . 

VII.  80GGE8TI0BS  FOR  FURTHER  WORK 

Since  for  the  ease  of  a  turbine  noszle,  the  aost  interesting  proper^ 
is  the  exit  Telocity  at  the  exit  area  of  the  nozzle,  a  aeasured  Telocity 
distribution  along  the  nozzle  would  be  very  useful.  Due  to  the  fact  that 
condensation  occurs  (mly  in  the  siq)ersonic  pert  of  the  nozzle  velocity 
aeasureaents  azr  quite  difficult  for  the  very  saall  dimensions  of  the 
nozzle. 

Tenqterature  aeasureaents  could  be  done  either  by  optical  methods  or 
by  using  themo-coiqples,  but  even  for  this  problem  the  dimensions  of  the 
nozzles  xised  are  too  small. 

For  further  work  to  study  velocity,  tenqperature  and  pressure  distrib¬ 
utions  which  differ  from  the  simplified  one-dimensional  case  a  much  bigger 
nozzle  is  reconnended  and  this  would  require  a  closed  cycle  with  a  condenser 
after  the  nozzle,  followed  by  a  conqpressor.  In  this  case  the  minimal  tem¬ 
perature  at  the  low  pressure  side  would  give  trouble  in  the  condenser. 
Therefore,  the  whole  cycle  would  probably  demand  a  higher  piressure  level. 

It  might  be  possible  to  reduce  the  scatter  by  obtaining  a  more  exact 
shape  of  the  pressure  distribution,  especially  near  the  discontinuity. 

This  could  be  done  by  using  a  search  tube  instead  of  pressure  taps,  so 
that  the  pressure  points  could  be  made  as  closely  spaced  as  possible. 

Perhaps  direct  measiurements  of  actual  size  of  condensed  droplets,  by 
light  scattering  or  other  techniques,  could  lead  to  a  better  understanding 
of  the  phenomenon  of  drop  growth,  and  to  a  better  appraisal  of  the  theory. 


In  toy  case,  a  better  teaperature  control  thui  the  one  adopted  in  these 
esqperiaents  would  be  a  desirable  isqprovement . 
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C  JRTRAN  PROGRAM  FOR  EXPANSION  OF  AMMONIA  IN  NOZZLES 


MAIN  PROGRAM 
READ  DATA 
READ  1»ANO,BOLK,G 

1  FORMAT  (2E14,6»F10.0) 

READ2.A.TK,EOT 

2  FORMAT  (3F10.0) 

READ4,DELX»DELMA 

4  FORMAT(2FlO.O) 

READS.  DELM.DELTAU 

5  FORMAT  (2F10.0) 

READS.  DELR.XI 

6  FORMAT  (2F10.0) 

READ  667.AJKR 
READ  667. FACTOR 

667  FORMAT  (E14.6J 

READ  67.ASTAR.XSTAR 
67  FORMAT  (2F10.0) 

READ  7.M 

7  FORMAT  (110) 

DO  777  I=1«M 

777  READ  77.XT( I ) .PADT( I ) 

77  FORMAT  (2F10.0) 

READ  7001.TOL(1) .T0L(2) .T0L(3) .T0L(4) .TOL ( 5 ) . TOL ( 6 ) ♦ TOL ( 7 ) 

7001  FORMAT  (7F10,0) 

OREAD  7002.YMIN(1 ) . YM IN ( 2 ) » YMI N ( 3 ) . YMJ N ( ^ ) ♦ YM I N ( 5 ) . YMIN ( 6 ) t YMI N ( 7 ) 

7002  FORMAT  (7F10.0) 

READ  7003.L.DXOUT 

7003  FORMAT  (I10.F?0.0) 

READ  7004.MARk(2) 

7004  FORMAT  (110) 

READ  8. BETA 

8  FORMAT  (FIO.O) 

33  REA03. PZERO. TZERO 

3  FORMAT  (2F10.0) 

READ  81. NUMBER 

81  FORMAT  (ID 
C  PART  1 

PRINT  298. PZERO 

298  format  (IOh  PZERO  »  E20.8) 

PRINT  299. TZERO 

299  FORMAT  (lOH  TZERO  =  E20*8) 

PRINT  300 

3000FORMAT( IlOH  X  AC  M 

IP  V  T  ///) 

COMMON  BETA.M. I .XT.RAOT 

COMMON  TK.EOT.A.R.ANO.BOLK.G.DELR.XI .DELTAU.AMASS  .FACTOR 
DIMENSION  I(50).XT(50) .RADT(50) 

R=ANO*BOLK/A 

T=T2ERO 

XsO.O 

10  X=X+DELX 


CALL  N022LE(X.AC.DADX) 

CP»0.4212+0.0151*T/100.+0,0038051*(T/100.)*( T/100. ) 

CP*42700.*CP 

AK*CP/(CP-R) 

A1=<AK-1. ) /(Af  +1*) 

A2=(  AC/ASTAR)ir*(2.*Al ) 

IF  (X-XSTAR)  15»16»17 

16  AMACH=1.0 
60  T025 

15  AMACH=r,0 

200  YMACH=A2»AMACH**(2.*A1)-A1*AMACH**2  -2./(AK+1,) 

IF  (YMACH)  56t2S*58 
56  XA=AMACH 

AMACH=AMACH+DELMA 
GO  TO  200 
58  XB=AMACH 

AMACH=0.5*{XA+XB) 

GO  TO  20 

17  AMACH=1.0 

2000  YMACH=A2*AMACH**(2.*A1 )-A1»AMACH**2  -2./(AK+le) 

IF  (YMACH)  158.25,156 
156  XA=AMACH 

AMACH=AMACH+DELMA 
GO  TO  2000 
158  XB*AMACH 

AMACH=0.5*(XA+XB) 

20  YMACH=A2*AMACH**(2.*A1 )-Al*AMACH«»2  -2,/(AK+l,) 
DYMACH=2.*A1*(A2*AMACH**( (AK-3.)/(AK+l,) J-AMACH) 
AMAC=AMACH-YMACH/DYMACH 
IF{ABSF(AMAC-AMACh)-DELM)  21.21.22 

22  AMACH.AMAC 
GO  T020 

21  AMACH*AMAC 

25  DEN0M  =  1,+0.5*(AI^-1.  )*AMaCH**2 
T«T2ER0/DEN0M 

P*P2ER0/DEN0M**(AX/ I AK-1. ) ) 

V=AMACH*SaRTF ( AX*G*R*T ) 

PRINT  69.X,AC.AMACH.P.V.T 
69  FORMAT  (6E20.8) 

0  PLG»9,584586-1648.6068/T-1,638646*T#0.01 
1  +2.403267*T*T#0.00001-1.168708*T»T*T*0.00000001 
PSIN=EXPF(2.302585*PLG) 

PSIN»1000,*PSIN 
IF  (P-PSIN)  If  ,10.205 
C  PART  2 

205  RHO  =  235.  +  0.654*(TX-T)  +  62 .77* ( Tk-T ) **0 . 333333 
RHOG  *  235.  +  0.654*(TK-T)  -  62 . 77* ( TK-T ) **0. 333333 
RHO  =  0.001*RHO 
RHOG  =  0.001*RHOG 

SI6IN  s  EOT*(TX-T)/'<A/RhO)**0. 666667 

HFG  =  32.938*SQRTF(406.16  -  T)  -  0. 5890*  ( 406 . 16  -  T) 

HFG  =  42700. *HFO 

UFG  =  HFG  “  PSIN*(1./RH0G  -  l./RHO) 

CP  »  0.4212  +  0.()151*T/100.  0.0038051*T*T/10000. 

CP  =  CP*42700. 
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AK  =  CP/(CP-R) 

AM  =  A/Af^O 

AV  *  AM/RHO 

RM  =  0,6203bl*AV**0. 333333 

Cl  =  0.797885*(PSIN/T)**2  »AV/BOLK  «SQRTF ( G/AM ) /BOLK 

C2  =  4,18879/(B0LK*T) 

AY  =  0.5*RHO*R*T»LO6F(P/PSIN) 

SIG  =  SIGIN  -  DELR»R.M*AY 

AJ  =  C1*(P/PSIN)**2  *  SQRTF{6IG)*EXPF(-C2*c>IG**3/AY**2) 
AJ=AJ#FACTOR 

IF  (AJ  -  AJKR)  100»  100*  ADQ 
100  GO  TO  10 

400  N  =  7 

PRINT  310 

310  F0RMAT{22H  CONDENSATION  OCCURS) 

PRINT  311»AJ 

311  FOaMAT(7H  AJ*  E20,8) 

PRINT  320 

3200r0RMAT{ 102H  X  Z1  Z  S 

1  MU  P  V  T/) 

C  AMASS-C0MPU7ATI0N 

PSIMAX=(2./(AK+1.  )  )**(1./<AK-1,)  )*SQRTF(A»v/(AK+l,  )  ) 
>.MASS=AS7AR#PS  IMAX*PZFR0*SQRTF  {2.*G/  (  R*TZER0)  ) 

Yd)  =  0.0 
Y{2)  =  0.0 
Y(3)  =  0.0 
Y{4)  =  C.C 
Y(5)  =  P 
Y{6)  =  V 
“(7)  =  T 
n  =  OELX 
XOUTd)  =  X  ' 

00  401  K=2»L 

4  ;1  XOUT(K)  =  X0UT(K  -  1)  +  DXOUT 
MARKd)  =  1 
MARK (3)  =  3 
MARK(4)  =  0 
MARK(5)  =  1 

DIMENSION  Y(7).  YMIN(7)»  TOL(7)*  XOUT(200)f  MARK(5) 

CALL  RUNGEINf  X,  Y,  TOL,  YMlNt  H*  XOUT,  MARK) 

GO  TO  (33*30) I  NUMBER 
30  CALL  EXIT 
END 


SUBROUTINE  DiFFEQ 
SUBROUTINE  DIFFEQ  (N.X»Y*0Y) 

DIMENSION  Y(7) »DY{7) 

COMMON  UETA.M, I ,XT»RADT 

COMMON  tk.eot.a»r.ano»bolk*g*oelr*xi,deltau»amass 

DIMENSION  I  (50  )  »XT(50) .RADTISO) 

P*Y(5, 

V»Y(6) 

T»Y( 7) 


»FACT0R 


A.i> 


CALL  NCZZLE(Xf ACtDADXJ 

2000PL&=9. 584586-1648. 6068/T-l. 63 8646*T*0» 01 

1  +2. 403267*T*T*0. 00001-1. 168708*T*T»T*0. 00000001 
PSIN  =  EXPF(2.5025851*PLG) 

PSIN  =  1000.*PSIN 

205  RhO  =  235.  +  0.654*(TK-T)  +  62. 77*( TK-T )**0. 333333 
RHO  =  0.001*RHO 
RH0G  =  PSIN/(R*T  ) 

SIGIN  =  E0T*( TK-T)/<A/RHO)**0. 666667 

HFG  =  32.938»SORTF(406.16  -  T)  -  0.5890*1406.16  -  Tj 
HFG  =  42700. *HF6 

UFG  =  HFG  -  PSIN*(1 ./RHOG  -  l./RHO) 

CP  =  0,4212  +  0.0151*7/100.  +  0.0038051*7*7/10000. 

CP  =  CP*42700. 

AK  =  CP/(CP-R) 

IF  (P/PSIN-1.0)  300. 300.301 

300  AJ=0.0 
RZERO=0.0 
GO  70  60 

301  AM  =  A/ANO 
AV  =  AM/RHO 

RM  =  0.620351*AV**0. 333333 

Cl  =  0.797885*(PSIN/7)*»2  *AV/BOLK  *SQR7F(G/AM) /80LK 
C2  =  4.18879/(80LK*7) 

AY  =  0.5*'RHO#R*7*LOGF(P/PSIN) 

SI6  =  SIGIN  -  OELR*RM*AY 

AJ  =  C1*(P/PSIN)**2  »  SQR7F{SIG)*EXPF{-C2*SIG**3/AY**2) 

AJ=AJ*FAC70R 

RS7AR  =  SIG/A' 

RZERO  =  1.3*RS7AR 
60  ALAM  =  HFG/{CP*7) 

AMACH  =  V/SQR7F(AK*G*R*7) 

IF(Y(3) )70.70»71 

70  RDROP*RZERO 
GO  70  122 

71  RDROP  =  3.*AMAS3«Y(4)/(RH0*Y(3)  ' 

122  CALL  7EMP(P. 7.70. RDROP) 

DRDX  =  1.5*P*R*(7D-7)/(V*UFG*Rh0*SQR7F(6.28319*R47  /G)) 

DY(1)  =  25.13274*AJ*AC 

DY(2)  =  Y(1)*DRDX  +  DY(1)*RZER0 

DY(3)  =  Y(2>*DRDX  +  12.56637*AJ*AC*RZER0**2 

DY(4)  =  RhO*Y( 3)*DRDX/AMASC  +  RH0*AJ*AC*4 ,1887904*RZERO**3/AMASS 
DY(5)  =  Y(5)*AK*AMACH**2/(AMaCH**2  -  1.)  *  U ALaM-1 , / ( 1 ,-Y ( 4 ) ) ) 

1  ♦  0Y(4)  -  DADX/AC) 

DY{6)  =  -Y(6)*DY(5)/(A<*AMACH**2*Y(5)  ) 

DY(7)  =  Y(7)*(AK-1. )/AK  *  DY(5)/Y(5)  +  ALAM*DY(4)  *Y(7) 

RE7URN 

END 


SUBR0U7INE  NOZZLE  FOR  SQUARED  NOZZLE’ 
SUBROU7INE  NOZZLE  (X.AC.DADX) 

COMMON  SE7A 


IF  (X-2.540)  11»12»12 

11  ROOT=SQRTF(5.f 8*X-X*X) 

AC*2. 770-1. 01b*ROOT 
DADX=-1.016*(2.54-X) /root 
GO  TO  300 

12  ASIN=SINF(BETA) 

IF  (X-2.540-2.540*ASIN}  11.21.21 
21  ACOS=COSF<BETA) 

0AC=0. 1885+2, 581»(l.-ACOS)  + 

1  1.016*(X-(2. 540+2. 540*ASINJ )*ASIN/AC0S 

DADX=1.016*ASIN/ACOS 
300  RETURN 
END 


C  SUBROUTINE  NOZZLE  FOR  CORRECTED  NOZZLES 

SUBROUTINE  NOZZLE ( X .AC .D^OX )  — — 

COMMON  BETA 
AIC«1.31 
ASTAR=et:TA 
A1«(AK-1.)/AK 
A2s{2,*AK-l, )/AK 
A3e(AK+l.)/{2.*AK,-2,) 

A4*(3.-AK)/(2.*AK-2,) 

CALL  TABLE(X.RAD.DRAD) 

AMACH=S0RTF( 2 .♦( 1 , /RAD**A1-1 , ) / ( AK-1 . ) ) 
DMACH—ORAD/  (  AK*AMACH*RAD**A2  ) 

B*(2,+(AK-1, )*AMACH*AMACH)/{AK+1. ) 
AC«ASTAR/AMACH*B#*A3 

0ADX*ASTAR4DMACH*(B**A4-B#*A3/(AMACH*AMACH) ) 

RETURN 

END 


C  SUBROUTINE  TABLE 

SUBROUTINE  IABlE  U.RAO.ORAD) 
COMMON  BETA.M. I .XT.RADT 
DIMENSION  I ( 50 ) .XT(50) .RA0T(50) 
1  =  1 

70  IF  (X-XT( I )  )  71.72.73 

72  RAD=RADT(  I  ) 

X1»XT( ■ ) 

X2  =  XT( I+l  ) 

RAD2  =  RADT(  I  +  U 
DRAD=(RAD2-RAD)/(X2-X1  ) 

GO  TO  69 

73  X1«XT(I) 

RAD1  =  RADT ( I  ) 

1  =  1  +  1 

GO  TO  70 

71  X2«XT(n 
RAD2=RADT( I ) 

DRAD=(RAD2-RAD1 )/(X2-Xl) 


RAD=RAD1+DRAD» (X“X1 ) 
69  CONTINUE 
RETURN 
END 


C  SUBROUTINE  TEMP 

SUBROUTINE  TE^P ( P » T t TD *RDR0P ) 

COMMON  BETA»M» I tXTtRADT 

COMMON  TK,EOT,AtR.ANOtBOLK.GtDEL.R.XI  ♦DELTAU.AMASS 
DIMENSION  I (50) »XT(50) .RADTtSO) 

DELTA=100.*DELTAU 

T1=T 

T2  =  T 

TAU=1. 

2000PLG=9 • 584586-1 648. 6068/T-l. 63 86A6*T*0. 01 

1  +2. 403267*T*T»0. 0000 1-1. 168708*T*T*T*0. 00000001 
PSIN  =  EXPF(2.302585j»PLG) 

PSIN  =  1000.*PSIN 

205  RHO  =  235.  +  0.6S4*(TK-T)  +  62 .77« ( TK-T J **0.333333 
RHO  =  0.001*RHO 
RHOG=PSIN/{R*T  ) 

SI61N  =  EOT*(TK-T)/ (A/RHO)**0. 666667 

HF6  =  32.938*SQRTF(406.16  -  T)  -  C.5890*{406. 16  -  T 

HFG  =  4270C.*HF'6 

UFG  =  HFG  -  PSIN*(1./RH0G  -  l./RHO) 

301  AM  =  A/ANO 
AV  =  AM/RHO 

RM  =  0.620351*/:v**0. 333333 
SIG=SIGIN/(1.+DELR*RM/RDR0P) 

122  B1  =  0.666667*UFG*XI/(R*T) 

B2  =  2.'^SIG/(RH0*R♦T*RDR0P) 

20  PI=(1.  .TAU-1. >/B1)*SQRTF(TAU)/EXPF(B2/TAU) 
PSINDR=PI*P 
CALL  SAT(PSINDR.TDR) 

T  =  TDR 

IF(A8SF(T2-TDR)-DELTA)  2*2.1 

1  r2=TDR 
TAU=TDR/T1 
GO  T0200 

2  TD=TDR 
T  =  T1 
RETURN 
END 


C  SUBROUTINE  SAT  ‘ 

SUBROUTINE  SAT (PS  I  NOR  .  TOR ) 

T=20.0 

DELT=50.0 

DELTOR=0.0001 

PSIN=PSINDR/1000. 

PLGDR=0.434294*LOGF(PSIN) 


A-7 


2000PL6=9. 584586-1 648. 6068/T-l, 638646*1*0. 01 

1  *2. 403267* T*T*0.0C001-1.168708*T*T*T*0. 00000001 
IF  (PL6-PLGDR)  10.22.30 
10  A=T 

T=T+DELT 
GO  TO  200 
30  B=T 

300  T*0.5*(A+BJ 

12000PLG=9. 584586-1 648. 6068/T-I. 638646*7*0.01 

1  +2. 403267* T*T*0. 00001-1. 16S708*T*T*T*0. 00000001 
0DPLG=164a.6058/(T*T ) -0.0 16 38646 
1  +2. 4032 67*0.00002*1-1, 168708*0. 00000003*T*T 

TI*T-{PLG-PLGDR)/DPLG 
IF(A6SF(TI-T)-DELTDR)  20»20.21 

21  T=TI 

GO  TO  1200 
20  T=T1 

22  TDR=T 
RETURN 
END 


5v^  \jr  MrrAKAIUS 


FIG  2a 

TEMPERATURE  BATH  AND  2-DIM- NOZZLE- 


FIG.  2b  SCHEMATIC  OF  STAGNATION  TANK  AND  AXI-SYMMETRIC 
NOZZLE  PROFILE 


FIG.  3  2- DIMENSIONAL  NOZZLE  PROFILE 
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FIG.  5  TRANSDUCER  CALIBRATION,  ILL’ STRATING  LINEAR  CHARACTERISTIC 
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FIG.  7  PRESSURE  D  ISTRIEiJT  I  ON  -  AX  I  -  SYMMETRIC  NOZZLE 


EXPERIMENTAL 


FIG.  8  PRESSURE  DISTRIBUTIONS,  CONTINUED,  2-DIMENSIONAL  NOZZLE  (I*  INCLUDED 
ANGLE) 
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PRESSURE  DISTR({3UT  ('.MS,  CONT«NUED,  2-DIMENSIONAL  NOZZLE  (r  INCLUDED 
ANGLE) 


FIG.  10  PRESSURE  DISTRIBUTIONS,  CONTINUED,  2-DIMENSIONAL  NOZZLE  0®  INCLUDED 
ANGLE) 


FIG.  It  PRKSSURE  DISTRIRUTIONS,  CONTINUED,  2-DIMENSIONAL  NOZZLE  (I*  INCLUDED 
ANGLE) 
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FIG.  12  PRESSURE  DISTRIBUTIONS,  CONTINUED,  2 -DIMENSIONAL  NOZZLE  (I®  INCLUDED 
ANGLE ) 
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FIG.  16  EFFECT  OF  Pq  ON  NON  -  CONDENSING  PRESSURE  DISTRIBUTION 
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FIG.  16  EFFECT  OF  Pq  ON  NON  -  CONDENSING  PRESSURE  DISTRIBUTION 


FIG.  17  PRESSURE  DISTRIBUTIONS  -  2  -  DIMENSIONAL  NOZZLE  (3®  INCLUDED  ANGLE) 
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FIG  18  ILLUSTRATES  UNSTEADINESS  NEAR 
NOZZLE  THROAT  AT  HIGH  R's 
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FIG.  20  THEORETICAL  PRESSURE  DISTRIBUTION  FOR  GEO 


FIG.  2IA  EFFECT  OF  CORRECTION  OF  SURFACE  TENSION  ON  THEORETICAL  PRESSURE 
DISTRIBUTION 
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FIG.  2IB  EFFECT  OF  VARIATION  IN  NUCLEATION  RATE  I  ON  THEORETICAL  PRESSURE 
DISTRIBUTION 
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FIG.  22  INCIDENCE  OF  CONDENSATION  IN  PRESSURE  - 
TEMPERATURE  -  DIAGRAM. 


